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Abstract: Recent phylogeographical studies indicated that glacial oscillations played a key 
role on the phylogeographic pattern of extant species. As most studies have previously 
been carried out on heavily ice-covered regions, such as in European and North American 
regions, potential effects of climatic oscillations on species that are distributed on ice-free 
regions are less known. To address this, we investigated the phylogeographic pattern of an 
avian species endemic to South China, which was not glaciated during the Pleistocene 
glaciations. By using 2142 bp mitochondrial DNA, we identified 89 haplotypes defined by 
39 polymorphic sites. A combination of high haplotype diversity (0.786-1.00) and low 
nucleotide diversity (0.00132-0.00252) was detected among geographic populations. 
Explicit genetic divergence was observed between S. s. semitorques and S. s. cinereicapillus 
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but not detected among geographic populations of S. s. semitorques. Divergence time of the 
two subspecies was dated back to 87 Kyr which is congruent with the interglacial MIS 5. A 
weak phylogeographic structure due to strong gene flow among geographic populations 
was identified in this species, suggesting complex topography of South China has not 
formed barriers for this species. 

Keywords: Collared Finchbill; East Asia; unglaciated region; phylogeography 



1. Introduction 

It is widely accepted that glacial oscillations during Pleistocene shaped the geographic distribution 
and genetic pattern of many plant and animal species on previously glaciated regions [1-3]. For species 
distributed in Europe and North America, their phylogeographic studies have shown congruent patterns. 
During glacial advance, these species receded to southern refugia. As ice sheets retracted, they 
experienced rapid expansions and recolonized northern latitudes after the Last Glacial Maximum 
(LGM) [3-5]. The repeated range expansions and retractions have considerably influenced their 
genetic diversity and population differentiations. 

In contrast, the influence of recent climatic oscillations due to glacial cycles on the phylogeographic 
patterns of species in unglaciated regions, especially in East Asia, remain poorly understood [3,6,7]. 
Recently, comparative phylogeographic research in an unglaciated North American region has shown 
congruent isolated and differentiation patterns with those of species distributed on previously glaciated 
areas, despite the fact that in some cases they may be older (Pliocene) [8]. 

East Asia harbors one of the most diverse temperate fauna as it was never heavily covered by 
ice-sheets [9,10]. During Pleistocene, East Asia was a mosaic of mountains lower than 2000 m and 
characterized by a relatively mild Pleistocene climate [11-14]. The climatic stability and 
heterogeneous topology of this region potentially hosted microclimatic zones capable of supporting a 
variety of habitats over time [15]. Three previous avian studies of this region showed different 
evolutionary history compared with Europe and North America [14,16,17]. Thus, glacial refugia may 
have been available for East Asian species throughout their entire ranges, instead of being limited to 
several regions of Europe and North America [14]. Accordingly, we assume that the milder 
Pleistocene climate may have evoked less stress for East Asian species than European and North 
American species, and possibly these species have similar patterns. 

To address this, we used two fragments of mitochondrial DNA to study the phylogeographical 
structure of a passerine species, the Collared Finchbill (Spizixos semitorques). The Collared Finchbill is 
endemic to south China, with two described subspecies: S. s. semitorques and S. s. cinereicapillus. 
S. s. semitorques cover a large area of South Mainland China, while S. s. cinereicapillus is only found 
in Taiwan. The Collared Finchbill distributes mainly in shrubbery, bamboo forests, and grassy areas on 
hilly and flat terrain below 1000 m, and its distribution is roughly between 20° N and 30° N in latitude, 
and 105° E and 120° E in longitude [18]. Understanding the genetic variation among populations and 
the past demographic processes sculpting the current population structure of S. semitorques might enrich 
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our understandings about the role of climatic changes during the Pleistocene glaciations on the 
phylogeographical formation of birds endemic to South China. 

2. Results and Discussion 

2.1. Results 

A combined length of 2142 bp mitochondrial genome, including Cytb (977 bp) and partial control 
region and tRNA-Glu (1165 bp) was obtained. Sequences aligned unambiguously have no gaps and 
there were no nonsense codons in protein-coding regions. Sixty- nine variable sites define 89 haplotypes. 
Most haplotypes are singletons, while there are only a few common haplotypes (10 see Attachment 1) 
shared by several geographical populations (Figure 1). The haplotype network depicts a compact 
topology and indicates no phylogeographic structure. 

Figure 1. Haplotypes median-joining network for Collared Finchbill generated by 
NETWORK 4.6. 




S. s. semitorques s - * cinereicapillus 



Neutrality test was conducted on geographic populations which contained at least five individuals. 
Tajima's D values detected negative values for all populations, but none are significant (Table 1). Fu's 
Fs statistics revealed significant negative values for SXFP (0.00) and GSWX (<0.01) populations. Fu 
and Li's D tests showed negative values for SXFP (<0.05), HNYJ (<0.01) and GZGD (<0.01) 
populations. No population showed significant negative value in all three neutrality tests. 
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Table 1. Neutrality test for geographic populations. 



S. s. semitorques S. s. cinereicapillus 



Statistics 


GSWX 


SXFP 


SCGY 


SCYA 


SNJ 


JXJGS 


HNYJ 


GZGD 


AHJX 


FJ 


TAIPEI 


Tajima's D 


-0.933 


-1.633 


-0.109 


-0.7064 


-0.996 


-0.694 


-0.591 


-0.691 


-0.986 


-1.199 


-0.894 


P 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


Fu's Fs 


-4.843 


-12.75 


-1.283 


-3.336 


-2.752 


-1.908 


-5.719 


-6.401 


-2.764 


-1.554 


-1.633 


P 


<0.01 


0.00 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


NS 


Fu and Li's D 


-1.313 


-2.803 


-0.109 


-0.640 


-1.270 


-0.798 


-0.765 


-1.360 


-1.370 


-1.199 


-0.894 


P 


NS 


<0.05 


NS 


NS 


NS 


NS 


<0.01 


<0.01 


NS 


NS 


NS 



Note: Three neutrality test parameters and statistical p values were calculated by DnaSP 5.0. 



AMOVA analysis show explicit divergence between two subspecies of Collared Finchbill (Table 2). 
We attempted to identify possible intraspecific genetic structure by defining several different ways of 
groupings in S. s. semitorques. However, AMOVA could not find a pattern with significant variation 
within any grouping. This indicates that there is no substantial genetic subdivision across subspecies on 
the mainland of China. Most of the variation is accounted for within-population variation, as could be 
expected given the large number of singletons. There is no correlation between F st genetic distance and 
geographic distance (Mantel test, Z = 1 18.0163, r = -0.1863, p = 0.9475) in S. s. semitorques (Figure 2). 



Table 2. AMOVA analyses of Collared Finchbill. 



Group 


Among 
groups 

(AG) 


Among 
populations 

within 
groups (AP) 


Within 
populations 

(WP) 


Percentage of variation 




FCT 


FSC 


FST 




2 Groups S. s. semitorques (mainland 
sites): S. s. cinereicapillus 
(TAIPEI) 


0.5043 


0.0098 


0.5092 


Among groups:50.44 
Among populations within 

group: 0.49 
Within populations: 49.08 


6 Groups East(AHJX, JSNJ): 
West(SCYA, SCGY): Southeast(FJ): 
Southwest(GZGD, GXJZS): Central 
(JXJGS, HNYJ, HNWG, SNJ): 
Northwest(GSWX, SXFP) 


0.0029 


0.0067 


0.0096 


Among groups: 0.29 
Among populations within 

group: 0.67 
Within populations: 99.04 



Note: For two groups, S. s. semitorques includes populations on mainland and S. s. cinereicapillus 
includes the population in Taiwan. 
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Figure 2. Data plot of the genetic distance (F st ) vs. geographical distance (Log geographical 
distance in km) of the Mainland Spizixos semitorques populations, showing the RMA 
regression line (R = 0.0167). 

030 -r- 




Geographic Distance (Log) 



MDIV gave an average of 0.44, 0.311 and 20.39 for TMRCA, t pop and 9 value respectively. 
Divergence time between two subspecies of Collared Finchbills was estimated to be about 87 Kyr. 
Meanwhile LAM ARC estimated 9 value for eight populations of S. s. semitorques. HNYJ population 
has a 9 value of 0.00806 which was larger than any population around it. ANJX population got 
relatively lower 9 value of 0.00249 and the remaining populations' results range from 0.00333 to 
0.00524. LAMARC also detected extensive gene flows among geographic populations of 
S. s. semitorques (Figure 3). 

Figure 3. Effective population size and gene flow among populations estimated 
by LAMARC. 
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2.2. Discussion 

Both local populations and the whole species of Collared Finchbill are characterized by high 
haplotype diversity and low nucleotide diversity. MtDNA markers detect an indistinctive genetic 
structure. These findings imply that the present populations of the S. semitorques seemed not to be 
affected by LGM and not segregated by topological complexity of South China. The results provide 
evidence of undiversified population structure. 

Genetic variation is likely homogeneously distributed among populations on mainland China. The 
result suggested an undistinguished genetic structure of S. s. semitorques haplotypes bear strong signals 
of: an irregular haplotypes network, high haplotype diversity and low nucleotide diversity, and no 
significant negative value of Tajima' D detected for all geographic populations. 

Results of AMOVA also show that there is no population subdivision on a regional scale, and only a 
small amount of variation could be attributed to groups (-1.5%). The majority of the variations (99% 
variation is found within populations, Table 2) were found within S. s. semitorques populations, 
possibly as a result of strong gene flows among geographic populations. Mantel Test analyses do not 
indicate the pattern of isolation by distance among all sampling sites of S. s. semitorques in South China. 
This indicates there is no geographic barrier to segregate S. s. semitorques. 

A pattern of high haplotype diversity and low nucleotide diversity is noted, which could be attributed 
to population expansion [19]. Given the harsh conditions during the Quaternary, there might have been 
numerous bottleneck events decreasing population size and genetic variability until such time as 
conditions became more favorable. 

Demographic History 

Although European and Asian continents cover similar latitude of Northern Hemisphere, they 
experienced different glacial histories [2,3,20]. Some compared phylogeographical researches 
indicated that Asian and European populations of the same species have obviously different genetic 
differentiation pattern [21-23]. This could reflect a different mechanism of species-specific responses 
to the glaciations. For example, many European species experienced post-glacial expansion due to the 
glaciations in the Quaternary, and expansion time roughly dated back to the retreat of the Last 
Maximun Glaciation (LGM) [2,3,20]. However, for many species distributes in Asia, their 
phylogeographic structures seemed not to be substantially affected by the LGM. Some extant species 
kept a stable population size during that time, but did not obviously expand with ice retreat after the 
LGM [14]. Three previous studies revealed intraspecific divergence and demographic history of some 
avian species in south China as a result of the Pleistocene climate change [14,16,17]. Mitochondrial 
data of Leucodioptron canorum canorum indicated the populations kept growing since penultimate ice 
age. This means complex topology of South China provided suitable habitats for L. c. canorum. 
Another example is Alcippe morrisonia, which widely distributes in South China and its ancestor dated 
back to 8 Myr. Geographical populations of A. morrisonia also kept stable or grew slightly during the 
LGM. Bambusicola thoracica thoracica showed a geographical structure among populations with little 
gene flow. 
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Unlike Alcippe morrisonia (Aves: Timaliidae), which is widespread in southern China and often 
flocks in the medium and understory of tropical rain forests and subtropical broadleaf evergreen 
forests [24,25], S. semitorques usually inhabits the canopy and edge of shrub land and has a larger body 
size [24] which might enable S. semitorques to have vagility and adaptability superior to A. morrisonia. 
The phylogeographical structure of A. morrisonia showed a deep lineage divergence while 
S. semitorques only have a shallow differentiation among geographical populations. Palynological and 
palaeoclimatic data suggest that the vegetation in East Asia in 35 Kyr ago similar to observed today 
might be ascribed to high precipitation, especially at middle or low latitudes [26-28]. The present results 
imply that during the cyclical oscillations of ice age, S. semitorques might stay relatively stable in 
certain areas of southern China. 

East Asia has experienced several cyclical changes due to the sea level rising and falling. These sea 
level fluctuations started in the Early Miocene (24 Myr BP), and became more frequent after 10 Myr, 
and particularly during the Quaternary. Such changes of sea level led to occasional 
connection/disconnection between mainland and Taiwan. Divergence time of the two subspecies dated 
back to 87 Kyr BP while the TMRCA was estimated to 0.44 Myr BP. This result implied that the rising 
sea level did not effectively block the two subspecies for a long time until one of the warmest interglacial 
marine isotope stage 5 (MIS 5:130-74 Kyr) [29]. The sea level was even higher at that time than it 
presently is [30]; so it is probably that the Taiwan Strait virtually blocked gene flow between two 
subspecies of S. semitorques from then on. 

3. Materials and Methods 

3.1. Sample Collection, DNA Extraction, Amplification and Sequencing 

A total of 120 samples of Collared Finchbill were collected from 14 localities from 2003 to 2009, 
covering most of its distribution range in southern China (Table 3 and Figure 4). Liver/muscle samples 
were stored in 100% ethanol immediately after removal. Total genomic DNA was extracted from tissue 
samples using the DNeasy Blood and Tissue Kit (QIAGEN) following the manufacturer's instructions. 
A partial cytochrome b (Cytb) gene of 977 (base pairs, bp) was amplified with the primer pair OSCL1 
and OSCH2 [31]. The thermocycling program consisted of an initial denaturation at 94 °C for 5 min, 
followed by 40 cycles of 94 °C for 40 s, 53 °C for 40 s and 72 °C for 40 s, plus a final extension at 72 °C 
for 5 min. We used the primers L437 and HI 248 [32] to amplify a fragment of 1 165 bp from the central 
D-loop region to the end of the tRNA-Glu gene, including tRNA-Pro and ND6. The PCR thermal profile 
started with 94 °C for 3min followed by 30 cycles of: 94 °C for 30 s, 55 °C for 30 s, and 72 °C for 30 s. 
Sequencing using the HI 248 primer was problematic, so a full sequence was obtained using two forward 
primers (16042: 5'-GTCACCAACTCCCAAAGC-3' and 16520: 5 '-ACCCAAAGCAAAACATAAACC-3 ') 
under designed using Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA, USA). All 
individuals were sequenced in both directions (ABI377). The sequences were deposited to GenBank 
(cytb: JF509464-509583 and D-loop to tRNA-Glu: JF509584-509703). 
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Table 3. Sampling information for Collared Finchbill. 



Subspecies 


Localities 


Latitude 


Longitude 


Sample 

size 


Haplotype numbers 


Nucleotide 
diversity (n) 


Haplotype 
diversity (Hd) 




GSWX 


32 5637.74 


104 41 0.16 


9 


U 1 3/14/15/ 16/17/18/19/21 


0.00236 


1.000 




SXFP 


33°11'24.0" 


108°12'0.0" 


25 


J_/l 1/14/20/24/26/27/28/29/30/31/32/ 
33/34/35/36/38/49/50/5 1 


0.00231 


0.980 




SCGY 


32°25'44.51" 


105°54'4.61" 


5 


7/8/9/10/12 


0.00243 


1.000 




SCYA 


29°35'24.0" 


102°35'24.0" 


10 


14/35/5 2/56/57/5 8/59/60/62 


0.00252 


0.978 




SNJ 


31°24'35.6" 


110°33'18.5" 


11 


2/3/4/5/6/37/45/55/61 


0.00236 


0.964 


5. s. 


HNYJ 


28°55'53.0" 


112°17'39.7" 


10 


22/24/25/42/43/44/46/47/48 


0.00240 


0.978 


semitorques 


HNWG 


26°25'48.0" 


110°22'12.0" 


3 


15/23/24 


0.00187 


1.000 




GZGD 


26°35'4.96" 


107°14'3.46" 


10 


1/15/49/68/69/70/7 1/72/73/74 


0.00245 


1.000 




GXJZS 


24°39'59.36" 


104°52'17.8" 


2 


37/67 


0.00187 


1.000 




JXJGS 


26°32'38.44" 


114°8'50.45" 


12 


24/75/76/77/78/79/80/8 1/82/83/84 


0.00195 


0.985 




JSNJ 


32°1'48.00" 


118°27'36.0" 


4 


53/54/56/57 


0.00210 


1.000 




AHJX 


30°4'15.12" 


118°35'44.45" 


9 


15/63/64/65/66 


0.00132 


0.786 




FJ 


27°19'54.86" 


118°7'13.85" 


5 


15/39/40/41/70 


0.00205 


1.000 


S. s. 

cinereicapillus 


TAIPEI 


24°54'56.56" 


121°40'26.17" 


5 


85/86/87/88/89 


0.00196 


1.000 



Note: Haplotype numbers are also used in network analysis and underlines denote shared haplotype. 



Figure 4. Sampling sites of S. s. semitorques: GSWX SXNS SCGY SCYA SNJ GZGD 
HNYJ JSNJ AHJX GXJZS JXJGS HNWG FJ; S. s. cinereicapillus: TAIPEI. Circles' size is 
proportional to amount of each sampling site. 



95'0'0 - E 100"0'0-E WOVE 110'0'0 - E 115°0'0 - E 12CrO'0"E 
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3.2. Sequence Analysis 

Fourteen sampling sites were identified for S. s. semitorques and one location for 
S. s. cinereicapillus (showed in Table 1). Sequences were aligned and edited using BioEdit 7.0 [33] and 
refined manually. DnaSP 5.0 [34] was used to calculate the number of segregating sites, haplotypes, 
haplotype diversity (h) [35], nucleotide diversity [36]. Genetic distance was computed using 
Arlequin 3.11 [37]. Tajima's D [38] was used to examine the selective neutrality of mitochondrial 
fragments. Two additional neutrality tests, Fu and Li's D [39] and Fu's F s [40] were used to detect 
departures from the mutation-drift equilibrium. Fs value tends to be negative if there is a significant 
excess of rare haplotypes, thus is often taken as evidence of recent demographic expansions or 
population bottlenecks [40]. All three tests were implemented in DnaSP 5.0 [34]. 

We used the combined sequences to reconstruct median-joining (MJ) network by program 
Network4.6 [41]. Mantel tests for isolation by distance (IBD) were performed using IBDWS [42] 

MDIV [43] was employed to estimate divergence time and migration rate between two subspecies. 
Aligned sequence data from sample sites were used to estimate the parameters H (=2N e l), 
M (=N e m = number of migrants between populations per generation), and T (the divergence time 
between populations, where 1 time unit = N e generations). Each analysis used an infinite-site model 
and 3000,000 generations of Markov chain Monte Carlo run and the first 500,000 generations was 
discarded as burn-in. Priors of maximum of M and T were set to 5 and 10 respectively. The divergence 
time between subspecies was estimated using the Formula divergence time = T pop x (Theta/2 uk). 

Molecular clock is widely used for the liable of sequence information. But the constancy of clocks 
have been controversial, particularly if they clash with estimates taken from more traditional sources 
such as the fossil record. Songbirds (Passeriformes) have a very limited fossil record in spite of 
making up more than one-half of the world's 10,000 living species of birds. There has been no fossil 
reported for Pycnonotidae birds. Different species have different molecular clock due to life cycles and 
metabolic rates. So for a specific species, a useful clock should be close to calibrated relatives. We 
used conventional 2% [44-46] molecular clock as mutation rate and 1.7 years [47] as generation time 
for Collared Finchbill. 

We estimated effective population size (0) and migration rate (M) in the assumption of both 
populations and migration are at equilibrium with the MCMC chain approach implemented in the 
program LAMARC 2.1.3 [48]. All analyses used a GTR model of evolution suggested by Model Test 
3.7 [49]. Two simultaneous searches with heating scheme were run 10,000,000 steps with an initial 
burn-in of 1,000,000 steps. LAMARC assumed that a population has been growing at the same 
exponential rate for a long period of time. Therefore, the ability to detect population growth depends on 
the amount of time a population experienced growth. Positive values of g indicated a growing 
population, and negative values indicated a shrinking population from past to present. The maximum 
probability estimates (MPE) for population under growth is above zero on average. As g value tended to 
be inflated, 95% CI of g including zero suggested undistinguished deviation from constant population 
size. Analyses were repeated three times with different random number seeds to assess consistency. The 
population growth or decline was significant if the 95% confidence interval of g did not include zero. 
The posterior probability densities and ESSs for these estimated population parameters were calculated 
in Tracer 1.5 [50]. 
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4. Conclusions 

Our study showed undistinguished genetic structure among geographical populations of 

5. semitorques. Approximate random gene flow indicated there might be no distinctive geographic 
barriers for this species in south China during late Pleistocene, which differs from the other three species 
of previous studies. Shallow differentiation due to considerable gene flow and lack of barriers might lead 
to high haplotype diversity and relative low nucleotide diversity. The mtDNA data showed a short 
evolutionary history of S. semitorques and undistinguished effects of the LGM on its 
historical demography. 

Acknowledgements 

Many thanks to Olivier Fran§ois for inviting us to contribute a paper for consideration and possible 
publication in a special issue on "Advances in Molecular Ecology". We thank Ornithological Research 
Group colleagues for their help with sample collection. This study was supported by the National 
Science Fund for Distinguished Young Scientists (No. 30925008), Major International (Regional) Joint 
Research Project (No: 31010103901) and CAS INFO-115-D02, NSFC (30870270, 30670276) to FM 
Lei, and by a NSFC program grant (J0930004) to Bin Gao. 

References 

1. Avise, J. Phylogeography: The History and Formation of Species; Harvard University Press: 
London, UK, 2000; pp 145-208. 

2. Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 2000, 405, 907-913. 

3. Hewitt, G. Genetic consequences of climatic oscillations in the Quaternary. Phil. Trans. Roy. Soc. 
London B 2004, 359, 183-195 

4. Taberlet, P.; Fumagalli, L.; Wust-Saucy, A.G.; Cosson, J.F. Comparative phylogeography and 
postglacial colonization routes in Europe. Mol. Ecol. 1998, 7, 453-464. 

5. Lessa, E.P.; Cook, J. A.; Patton, J.L. Genetic footprints of demographic expansion in North 
America, but not Amazonia, during the Late Quaternary. Proc. Nat. Acad. Sci. USA 2003, 100, 
10331-10334. 

6. Zhang, Y.; Ge, S. Molecular evolution study in China: Progress and future promise. Phil. Trans. 
Roy. Soc. London B 2007, 362, 973-986. 

7. Qiu, Y.X.; Fu, C.X.; Comes, H.P. Plant molecular phylogeography in China and adjacent regions: 
Tracing the genetic imprints of Quaternary climate and environmental change in the world's most 
diverse temperate flora. Mol. Phylogenet. Evol. 2011, 59, 225-244. 

8. Soltis, D.E.; Morris, A.B.; McLachlan, J.S.; Manos, P.S.; Soltis, P.S. Comparative phylogeography 
of unglaciated eastern North America. Mol. Ecol. 2006, 15, 4261-4293. 

9. Shi, Y. Quaternary glaciation in China. Quatern. Sci. Rev. 1986, 5, 503-507. 

10. Liu, K. Quaternary history of the temperate forests of China. Quatern. Sci. Rev. 1988, 7, 1-20. 

11. Weaver, A.; Eby, M.; Fanning, A.; Wiebe, E. Simulated influence of carbon dioxide, orbital 
forcing and ice sheets on the climate of the last glacial maximum. Nature 1998, 394, 847-853. 



Int. J. Mol. Sci. 2011, 12 



3645 



12. Pinot, S.; Ramstein, G.; Harrison, S.; Prentice, I.; Guiot, J.; Stute, M.; Joussaume, S. Tropical 
paleoclimates at the last glacial maximum: comparison of Paleoclimate Modeling 
Intercomparison Project (PMIP) simulations and paleodata. Clim. Dynam. 1999, 15, 857-874. 

13. Ju, L.; Wang, H.; Jiang, D. Simulation of the last glacial maximum climate over East Asia with a 
regional climate model nested in a general circulation model. Palaeogeogr. Palaeoclimatol. 2007, 
248, 376-390. 

14. Li, S.; Yeung, C. Sailing through the Late Pleistocene: unusual historical demography of an East 
Asian endemic, the Chinese Hwamei (Leucodioptron canorum canorum), during the last glacial 
period. Mol. Ecol. 2009, 18, 622-633. 

15. Qian, H.; Ricklefs, R.E. Large-scale processes and the Asian bias in species diversity of temperate 
plants. Nature 2000, 407, 180-182. 

16. Song, G.; Qu, Y.; Yin, Z.; Li, S.; Liu, N.; Lei, F. Phylogeography of the Alcippe morrisonia (Aves: 
Timaliidae): Long population history beyond late Pleistocene glaciations. BMC Evol. Biol. 2009, 
9, 143. 

17. Huang, Z.; Liu, N.; Liang, W.; Zhang, Y.; Liao, X.; Ruan, L.; Yang, Z. Phylogeography of 
Chinese bamboo partridge, Bambusicola thoracica thoracica (Aves: Galliformes) in south China: 
Inference from mitochondrial DNA control-region sequences. Mol. Phylogenet. Evol. 2010, 56, 
273-280. 

18. MacKinnon, J.; Phillipps, K.; He, F. A Field Guide to the Birds of China; Oxford University Press: 
New York, NY, USA, 2000; pp 362-363. 

19. Grant, W.; Bowen, B. Shallow population histories in deep evolutionary lineages of marine fishes: 
insights from sardines and anchovies and lessons for conservation. J. Hered. 1998, 89, 415-427. 

20. Hewitt, G. Some genetic consequences of ice ages, and their role, in divergence and speciation. 
Biol. J. Linn. Soc. 1996, 58, 247-276. 

21. Kryukov, A.; Iwasa, M.; Kakizawa, R.; Suzuki, H.; Pinsker, W.; Haring, E. Synchronic east-west 
divergence in azure-winged magpies (Cyanopica cyanus) and magpies (Pica pica)*. J. Zool. Syst. 
Evol. Res. 2004, 42, 342-351. 

22. Zink, R.M.; Drovetski, S.V.; Questiau, S.; Fadeev, I.V.; Nesterov, E.V.; Westberg, M.C.; Rohwer, S. 
Recent evolutionary history of the bluethroat (Luscinia svecica) across Eurasia. Mol. Ecol. 2003, 
12, 3069-3075. 

23. Pavlova, A.; Zink, R.M.; Drovetski, S.V.; Red'kin, Y.; Rohwer, S. Phylogeo graphic patterns in 
Motacilla flava and Motacilla citreola: species limits and population history. Auk 2003, 120, 
744-758. 

24. Zou, F.; Chen, G. A study of understory bird communities in tropical mountain rain forest of 
Jianfengling, Hainan Island, China. Acta Ecol. Sinica 2004, 24, 510-516. 

25. Taiwanica, A.Z. Diet analysis of the gray-cheeked Fulvetta (Alcippe morrisonia) at Fushan 
Experimental Forest in Taiwan. Acta Zool. Taiwan 1998, 9, 59-66. 

26. Yuan, D.; Cheng, H.; Edwards, R.; Dykoski, C; Kelly, M.; Zhang, M.; Qing, J.; Lin, Y.; Wu, J. 
Timing, duration, and transitions of the last interglacial Asian monsoon. Science 2004, 304, 
575-578. 



Int. J. Mol. Sci. 2011, 12 



3646 



27. Kelly, M.; Edwards, R.; Cheng, H.; Yuan, D.; Cai, Y.; Zhang, M.; Lin, Y.; An, Z. High resolution 
characterization the Asian Monsoon between 146,000 and 99,000 years BP from Dongge Cave, 
China and global correlation of events surrounding Termination II. Palaeogeogr. Palaeoclimatol. 
2006, 236, 20-38. 

28. Yu, G.; Gui, F.; Shi, Y.; Zheng, Y. Late marine isotope stage 3 palaeoclimate for East Asia: A 
data-model comparison. Palaeogeogr. Palaeoclimatol. 2007, 250, 167-183. 

29. Zhao, B.; Wang, Z.; Chen, J.; Chen, Z. Marine sediment records and relative sea level change 
during late Pleistocene in the Changjiang delta area and adjacent continental shelf. Quatern. Int. 
2008, 186, 164-172. 

30. Hearty, P.J.; Hollin, J.T.; Neumann, A.C.; O'Leary, M.J.; McCulloch, M. Global sea-level 
fluctuations during the Last Interglaciation (MIS 5e). Quatern. Sci. Rev. 2007, 26, 2090-2112. 

31. Qu, Y.; Lei, F. Comparative phylogeography of two endemic birds of the Tibetan plateau, the 
white-rumped snow finch (Onychostruthus taczanowskii) and the Hume's ground tit 
(Pseudopodoces humilis). Mol. Phylogenet. Evol. 2009, 51, 312-326. 

32. Tarr, C. Primers for amplification and determination of mitochondrial control-region sequences in 
oscine passerines. Mol. Ecol. 1995, 4, 527-530. 

33. Hall, T. BioEdit: A User- Friendly Biological Sequence Alignment Editor and Analysis Program 
for Windows 95/98/NT. Nucleic Acids Symp. 1999, 41. 95-98. 

34. Librado, P.; Rozas, J. DnaSP version 5.0: A software for comprehensive analysis of DNA 
polymorphism data. Bioinformatics 2009, 25, 1451-1452. 

35. Nei, M. Molecular Evolutionary Genetics; Columbia Universiy Press: New York, NY, USA, 1987. 

36. Nei, M.; Tajima, F. DNA polymorphism detectable by restriction endonucleases. Genetics 1981, 
97, 145-163. 

37. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for 
population genetics data analysis. Evol. Bioinform. 2005, 1, 47-50. 

38. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. 
Genetics 1989, 123, 585-595. 

39. Fu, Y.; Li, W. Statistical tests of neutrality of mutations. Genetics 1993, 133, 693-709. 

40. Fu, Y.X. Statistical tests of neutrality of mutations against population growth, hitchhiking and 
background selection. Genetics 1997, 147, 915-925. 

41. Bandelt, H.; Forster, P.; Rohl, A. Median -joining networks for inferring intraspecific phylogenies. 
Mol. Biol. Evol. 1999, 16, 37-48. 

42. Jensen, J.; Bohonak, A.; Kelley, S. Isolation by distance, web service. BMC Genet. 2005, 6, 13. 

43. Nielsen, R.; Wakeley, J. Distinguishing migration from isolation: A Markov chain Monte Carlo 
approach. Genetics 2001, 158, 885-896. 

44. Weir, J.; Schluter, D. Calibrating the avian molecular clock. Mol. Ecol. 2008, 17, 2321-2328. 

45. Weir, J.T. Divergent timing and patterns of species accumulation in lowland and highland 
neotropical birds. Evolution 2006, 60, 842-855. 

46. Lovette, I.J. Mitochondrial dating and mixed support for the "2% rule" in birds. Auk 2004, 121, 
1-6. 



Int. J. Mol. Sci. 2011, 12 



3647 



47. Saether, B.E.; Lande, R.; Engen, S.; Weimerskirch, H.; Lillegard, M.; Altwegg, R.; Becker, P.H.; 
Bregnballe, T.; Brommer, J.E.; McCleery, R.H.; Merila, J.; Nyholm, E.; Rendell, W.; 
Tryjanowski, P. Generation time and temporal scaling of bird population dynamics. Nature 2005, 
436, 99-102. 

48. Kuhner, M. LAMARC 2.0: Maximum likelihood and Bayesian estimation of population 
parameters. Bioinformatics 2006, 22, 768-770. 

49. Posada, D.C.; Crandall, KA. MODELTEST: Testing the model of DNA substitution. 
Bioinformatics 1998, 14, 817-818. 

50. Rambaut, A.; Drummond, A.J. Tracer vl.5, 2007. Available online: http://beast.bio.ed.ac.uk/Tracer 
(accessed on 29 March 201 1). 



© 2011 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



